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ABSTRACT 


A numerical model incorporating heat transfer theory” with 
cylindrical geometry has been developed to investigate the sublimation 
rate of dry ice pellets in clear air and simulated cloud. A_ series of 
thirty two experiments were conducted in the University of Alberta FROST 
icing-wind tunnel, and nine aircraft drop experiments were performed to 


_ verify model predictions of the sublimation rate. 


The principal conclusions drawn from this study are: (a) The FROST 
icing wind tunnel is an appropriate facility for conducting sublimation 
rate experiments, and may also be useful for investigations into ice 
crystal production rates. (b) In spite of the use of several 
simplifying assumptions, the cylindrical model predicts the sublimation 
rates of dry ice pellets to within 10%, when compared with wind tunnel 
observations. (c) The model predicts the fall times of cylindrical 
pellets to within 10%-15%, when compared with the free-fall experiments. 
(d) Cloudy and saturated conditions, with high temperatures, enhance the 
sublimation rate of dry ice but, cloudy and saturated conditions, with 
cold temperatures do not make an appreciable effect on the sublimation 
rate of dry ice. (e) Numerical simulations of freely-falling pellets 
suggest pellets with diameters of between 0.5 and 0.6 cm ffor efficient 


use of dry ice when dropped from the -10 °C level. 
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CHAPTER 1 


INTRODUCTION 


The manipulation of weather to his advantage has been a dream of 
man for many centuries. In its infancy, this ambition was limited to 
Tocal need (rain for a patch of corn), and the methods were primitive 
(rain dances and lighting fires, Hess,1974). As knowledge of 
atmospheric processes increased, man's approach to this dream was 
directed toward scientific consideration of the possibilities of 


altering these processes. 


The modern era of weather modification was initiated by Schaefer 
(1946), when dry ice (solid carbon dioxide), was serendipitously used in 
his efforts to discover a material to nucleate freezing in supercooled 
clouds. This event led to his discovery that dry ice can produce a host 
of ice embryos when inserted into a cloud. He observed that the 
concentrations produced by a tiny bit of dry ice were high and the ice 
crystals were very small. The mechanism by which dry ice produces large 
amounts of ice crystals in its wake while falling through a supercooled 


cloud was found to be homogeneous nucleation from the vapor phase. 


The recent revival of the use of dry ice as an alternative cloud 
seeding agent, (Schaefer,1976; Holroyd et al.,1978; Silverman,1979; 
English and Marwitz,1981), has renewed interest in the sublimation rate 


of dry ice pellets. Knowledge of the sublimation rate, to accurately 
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predict the distance that dry ice pellets will fall, is essential if the 
dry ice is to be used efficiently in cloud seeding operations. With 
such information, an optimum pellet size may be selected for a given 
mission,( provided that the pellets can be produced to whatever 
specified size ). The sublimation rate is also an important parameter 
for numerical modellers to use in their studies simulating the effects 


of dry ice seeding on cumulus clouds (Kopp et al.,1979). 


Some previous studies have been conducted to determine the 
sublimation rate of dry ice pellets. An investigation, (Mee and Eadie, 
1963), was undertaken to find a means of increasing the efficiency of 
the dry ice cloud seeding process, and to develop field techniques and 
equipment for U.S Army personnel to dissipate whiteouts. Several 
experiments were conducted in a vertical wind-tunnel, to establish the 
time required for the complete sublimation of dry ice pellets of a given 
initial mass falling at their terminal velocities. From these 
measurements an empirical expression for the total sublimation of dry 


ice as a function of initial mass has been proposed. 


Another study, (Fukuta, et al., 1971) was conducted to investigate 
the effectiveness (the number of ice crystals generated per gram of dry 
ice sublimed) of dry ice. Empirical equations for the Sublimation rate 
of dry ice, terminal velocity and pellet surface temperature were 
obtained. Both of these studies have assumed the pellet geometry to be 


spherical, and they were limited to sublimation in clear air. 


Summer cumulus cloud seeding . experiments (English and Marwitz, 
1981) have been conducted as part of the Alberta Hail Project since 


1978. Three seeding treatments: dry ice pellets, silver iodide flares, 
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and placebos are applied to these clouds. The dry ice pellets used in 


these experiments have been obtained commercially in cylindrical form. 


A few questions that needed to be resolved were: a) Wouldn't a 
sublimation rate expression based on cylindrical geometry be more 
appropriate to use for these experiments? b) Does a cloudy environment 
affect the sublimation rate of the dry ice? A study attempting to 
answer these questions was initiated, and the results are presented in 


this thesis. 


The approach taken in this thesis is one which incorporates both 
theory and experiment in an attempt to examine the sublimation rate of 
cylindrical dry ice pellets both in clear air and simulated cloud. A 


description of the theory and the model will be presented in Chapter 2. 


Laboratory experiments carried out in the University of Alberta's 
FROST icing-wind tunnel are discussed in Chapter 3. Their purpose was 
to test the the theoretical model. Descriptions of the FROST tunnel 
capabilities, the experimental procedures followed, and the data 
analysis techniques used are presented. The experimental results and 
comparisons with the theoretical model predictions of changes in 


diameter are also presented and discussed. 


Several field experiments were conducted in which cylindrical dry 
ice pellets were dropped from an aircraft flying at a known altitude. 
These experiments are discussed in Chapter 4. The numer ical model 
developed in Chapter 2 was modified slightly to investigate the 
sublimation of cylindrical dry ice pellets falling at terminal velocity 
through the atmosphere. The purpose of these experiments was to provide 


a check on the free-fall model predictions of the time taken to reach 
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the surface, and the final diameters the pellets achieve. A description 
of the free-fall model, and discussion of the results from these 


experiments are presented. 


Conclusions and recommendations for future work are presented in 


Chapter 5. 
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CHAPTER 2 


THE THEORETICAL MODEL 


2.1 Theoretical Model Description 


A theory is presented for the sublimation of fixed, cylindrical dry 
ice pellets in a moving, saturated (with respect to water) cloudy 
airstream. A heat flux balance (Lozowski et al. (1979)), at the 
surface of the cylindrical dry ice pellet is assumed. The heat balance 


equation can be expressed as: 
+ a = 
Q Q ue 1p Q. 7 Qy 0 (1) 


where Q. is the total rate of sensible heat transfer from the airstream 
to the pellet by conduction and convection (considered negative), on is 
the total rate of heat transfer to the pellet due to the sublimation of 
water vapor onto the dry _ ice (considered negative), Q. is the heat 
absorbed by the pellet due to the sublimation of the dry ice itself 
(considered positive), Qe is the net heat transfer due to radiation 
(considered negative), and Qy is the sensible and latent heat released 
due to the cooling and freezing of impacting water droplets (considered 


negative). 


End effects are ignored in the current version of the model. That 
is to say, sublimation is assumed to occur from the cylindrical surface 
only and not from the ends of the pellet. Furthermore, it is assumed 
that the sublimation takes place uniformly both around the circumference 


and along the length of the cylinder. Thus the pellet, in the model, 
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retains its cylindrical shape and a constant length as it sublimates. 
In order to strictly satisfy the cylindrical shape preservation 


assumption, the pellet would have to rotate about its longitudinal axis. 
The individual terms of equation 1 are described as follows: 


The total sensible heat transfer from the airstream to the 


cylindrical surface by conduction and convection is given bys: 
Q. = mDLh(t.-t,) (2) 


where D is the cylinder diameter, L its length, h is the _ heat 
transfer coefficient, and me and t, are, respectively, the pellet 
surface temperature and the airstream temperature. In the present work, 


the pellet surface temperature is taken to be -100 C (Fukuta, et al. 


1971). 


The heat transfer coefficient h is related to the Nusselt number 


Nu via the definitions 


k Nu 
pan a (3) 
Mi Nag 


where Kk. is the thermal conductivity of air at the ambient temperature, 


D is the cylinder diameter and Nu is the Nusselt number. 


The Nusselt number was calculated using an expression given by 


Zukauskas (1972) for smooth circular cylinders in a cross flow of air: 
Nu = 0.26 Re®® (4) 


This expression is valid over the range 10°< Re < 10°. which includes the 
Re range of practical interest. An alternate method, using the results 


from the wind tunnel experiments to determine a relation between the 
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Nusselt and Reynolds numbers was also explored. These results are 


presented in Appendix A. 


The total latent heat transfer due to sublimation of the water 


vapor onto the dry ice is given by (Lozowski et ale; 1979): 


+63 6] 
Qeg = ~mbth( EX) Fen (eg(t)e5(t,) (5) 

where Pr is the Prandt] number of air, SC is the Schmidt number for 
water veneer in air , e¢ is the ratio of the molecular weights of water 
vapor and dry air, le is the specific latent heat of sublimation of 
water at the surface temperature, P is the static pressure in the free 
stream, a is the specific heat capacity of dry air at constant 
pressure, ett) is the saturation vapor pressure of water vapor in the 
air, at the ambient temperature, and ett ) is the saturation vapor 


S 


pressure of water vapor at the surface temperature of the pellet. 


The saturation vapor pressure of water vapor at the surface 
temperature is so much less than the saturation vapor pressure at the 
ambient temperature, that it can be neglected and is assumed to be zero. 


The derivation of this expression is presented in Appendix B. 


The rate of consumption of latent heat required to give rise to a 


sublimation rate dm/dt (assumed negative) is given by: 
Q=-L, — (6) 


where m is the mass of the dry ice pellet, and L. is the specific 


latent heat of sublimation -of dry ice at the surface temperature. 
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The heat transfer due to radiation is given by: 


Q = -o(T! Sa Te e egdk OF) 


where o is the Stefan-—Boltzmann constant, and Ee, and ee are the 


emissivities of the surrounding air and the pellet surface. 


The heat transfer by the water droplets can be considered as a 
three stage process: the heat needed to cool the water droplets to 0°C, 
the heat needed to freeze the droplets at O0°C, and the heat needed to 
cool the frozen droplets to the pellet surface temperature. This can be 


expressed as: 
ae ee Senne oe 
Qy = (c ne Aaa O’C) + RY Le + Cc. R (0 C aay (8) 


where oe is the average specific heat capacity of water between t. and 
GC, RT is the mass flux of droplets collected per unit time, Le is the 
specific latent heat of freezing at 0°C, and C is the average specific 


heat capacity of ice between 0°C and toe 


The contributions to the heat balance due to the water droplets and 
radiation have been neglected in the model. The justification is found 


in Appendix C. 


Substituting expressions (2), (5), and (6) into equation (1) 
yields: 
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Substituting equation (3) into (9), expressing the mass as 
mD 2 Lo /4, where 0. is the dry ice pellet density (experimentally 
determined to be 1.4 gem’) and taking to be constant we obtain the 


following relationship: 


-63 
dD eal uy, Hes eee (t ) (10) 
dt P.L,D S a Sc ‘<= a’oa 


An analytic solution to this relationship exists when the ambient 
conditions (temperature and wind speed) can be held constant like in a 
wind tunnel. The solution is presented in Appendix 0D. However, a 
numerical technique was used to solve this’ relationship since the 
numerical model was also going to be used in simulating the free-fall 
experiments where conditions are not constant. Comparisons were made 
between the numerical solutions and the analytic solutions. The results 


were within 5%-10% of each other. 


The effect of relative humidity on the sublimation rate can be 
investigated by substituting the actual vapor pressure (the saturation 
vapor pressure multiplied by the relative humidity) for the saturation 
vapor pressure in equation 10. Measurements of the relative humidity 
within the test section both with the sprays off and on = were 
unfortunately never made. Consequently, the actual vapor pressure in 
the test section was not known. However, when the sprays were on the 
relative humidity was assumed to be 100%. Thus the One term was 


retained in cloudy conditions, but it was ignored in clear air. 
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2.2 Sensitivity Tests 


The most conmon way of numerically solving a differential equation 
is to find an approximate expression for the derivative appearing in the 
equation. These approximate expressions are formed using differences of 
the dependent variables over finite space and time intervals. This 
approach is called the finite difference method, and the resulting 
algebraic equation is called a finite difference scheme to that 


differential equation. 


A forward differencing finite difference scheme was used to _ solve 
equation 10. The model was run varying the time interval to check the 
numerical accuracy of the solution. The experimental conditions used 
were T= -10°C, V=l4ms7!, and time intervals of (Koy 2 le.0.5,. and 
0.25 sec. The results of these experiments are tabulated and appear _ in 


Table 1. 


si a , x 
i. i. = a ) a) : . 
i i a a mL ‘ 


f ' taatae | 
no \isupe tettneveYt: eo pnivioe remade Yenitt 1 vate poms — “— 
ent at gel “eeegs ov lfevl Se a2 407 nol 2eo7dxb! paige: et 
1o genre se?) | piel abs Arnot 2757 pieicins 70 mT seh sat: sesiit 
gist Lateretn! onto the ssct asbnts 220 Net aehaah 4 
‘onigtuebes ot? wre , bolsen in web ext cilS ang asada de < 
“gahs of ‘suptee sane witb, Stl = beidas q| ‘eiteatins 


Taitees Te 


de. Ras 


ue 
Bhi 1 a 


. a , 
awe CS fat ter iStetce a sheaeTti > aryl mT at T oceeceet) 1@ SUSHI Aye | 
wl) adnoeis Go [evisinh em’ aid” peice t+ ra eben mitt «Ot ; 
bem gnots: MIO” KI pon! 19%; <3 fo Pip ion, 607) 36 vos wats wv 
BAM. of wf “i S gt l,O? 320 2 hey ii) aim 3 ine ; “ame tay snit- 
‘al Serge bre ae ere gtrem! exe seen To 27 luca ey one 


4 ; ' , 


al 


At=0.25sec At=0.5sec At=1sec At=2sec At=Ssec At=10sec 

Time DIAMETER 
(sec) (cm) 

0 1.5 1.5 1.5 1.5 1.5 1.5 

60 1. 361454 1.361514 1.361562 1.361621 1.361766 1.361994 
120 1.217016 1.217147 1.207250.) 142107377 e217700: 1.216212 
180 1.065357 1.065560 1.065732 1.065948 1.066501 1.067381 
240 0.904547 0.904809 0.905050 .0.905377 0.906239 0.907621 
300 0.731349 0.731670 0.731997 0.732481 0.733798 0.735919 
360 0.539680 0.540097 0.540572 0.541323 0.543413 0.546778 
420 0. 314613 0. 315248 0.316063 0.317439 0.321308 0.327456 

Table 1 


predicted diameter after 7 minutes are within 4% of each other. 
time step increases from 0.25 sec to 10 sec 
marginally. 


difference scheme would be adequate, and all subsequent models were 


Sensitivity test results 


the error 


with this scheme, using a time step of 1 second. 


increases 


The results of these tests indicate that the differences in the 


These tests indicated that the forward differencing finite 
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CHAPTER 3 


THE WIND TUNNEL EXPERIMENTS 


3-1 FROST Tunnel Capabilities 


The icing wind tunnel is situated in the Department of Mechanical 
Engineering and was designed by M. Sroka (1972) and Prof. G. Lock. 
FROST is an acronym ffor Fundamental Research on Solidification and 


Thawing. 


A plan view of the icing wind tunnel is shown in Figure 1. The 
tunnel has an octagonal cross-sectional shape with a diameter of 0.46 m 
across the flat plates in the working section. The air is impelled 
through the tunnel by a 22 kW constant speed electric motor direct 
coupled to a fan. Velocity control is achieved by adjusting a set of 
vanes immediately downstream of the fan. The air is cooled by a 30 kW 
refrigerator system. The spray system consists of four pneumatic 
atomizing nozzles, capable of producing water droplets with diameters 


typically ranging between 10 Um to 100 UM, 


The icing wind tunnel is capable of producing air speeds up to 60 
ms~1 in dry air and 35 ms~'with the sprays on , liquid water contents up 
to 3 gm73, and temperatures down to -18 °C. The rapid accumulation of 
ice especially on the turning vanes (Gates, 1981) accounts for the 


difference in maximum air speed with the sprays on. 
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3.2 Exper imental Procedure 


In the experiments we measured the rate of dry ice sublimation’ in 
an airstream. A cylindrical dry ice pellet with a measured diameter and 
length was suspended by means of an insulated clamp in the center of the 
46 cm diameter test section of the tunnel. The cylinder projected at 
least 5 cm beyond the clamp and was oriented so that its longitudinal 
axis was perpendicular to the airstream. The airflow was then started 
and photographs were taken every 10 or 15 seconds looking along the 
cylinder axis, to document the size and shape change of the 
cross-section. An ASAHI Pentax camera equipped with a super 
multi-—coated MACRO Takumar 1:4 50mm lens was used to photograph all the 
experiments. The sublimation was followed to the point where the 
cylinder became so small that it broke free of the clamp and was swept 
down the tunnel. Because the cylinder was clamped in a_ fixed 
orientation, its cross-section did not remain perfectly circular as it 
sublimed. In fact, it eventually achieved a wedge shape, as can be seen 
by the analyzed cross-sections depicted in Figure 2. These results are 
from experiment 5, conducted at an initial temperature of 18 °C and = aan 
airspeed of 14 ms7! . The photographs were taken at 10 second 


intervals. The degree of enlargement is 1.45 the original size. 
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The analyzed cross-sections of a dry ice pellet. 
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In order to make a comparison with the theory , therefore, the 
cross-sectional area was measured and an equivalent circular diameter 
was calculated. Photographs were taken every minute in some of the 
earlier experiments, to check for changes in length. Sub sequent 
analysis revealed that the lengths did not change by more than 10%, from 
the original measurement of the length, during the entire duration of 
the experiment. (As an example; dur ing experiment 10 the diameter 
changed from 1.40 cm to 0.80 cm after 4 minutes of elapsed time, 


while the length changed from 7.80 cm to 7.15 cm after 4 minutes) 


In order to simulate the sublimation rate of dry ice pellets 
falling through a cloud, some experiments were conducted in the FROST 
tunnel with the sprays on. Figure 3 shows the pellet in the working 
section of the tunnel with the sprays on. The wake of the pellet can be 


clearly seen. 


The measurement of the liquid water content was accomplished by 
measuring the mass of ice accreted on a rotating cylinder (DeLorenzis, 


1980). 


The wind speed in the working section was determined by taking 
observations of the water filled U-tube manometer, which measures the 


pressure drop across the contraction, and using: 
Ap = 1/29,v" (ts) 


where AP is the dynamic pressure, V is the airspeed, and P, is the 
density of air evaluated at the temperature and the pressure (930 mb) of 


the working section. This formula assumes that the air is 
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incompressible which may introduce errors that are most likely small 


since the maximum wind speed was 31 ms7!. 


The temperature in the working section was measured before and 
after every experiment using a thermometer. Difficulties in maintaining 
a constant temperature for the duration of some experiments were due to 
the heat energy caused by the friction of the airstream on the tunnel 
walls. A temperature balance is achieved when the frictional heat 
energy production rate is in equilibrium with the heat transfer through 
the tunnel walls. This balance at times took longer than the time 


required to complete a series of experiments. 


A total of 32 experiments were conducted in the FROST tunnel 
varying the temperature and airspeed, and with the sprays on and off. 


The exper imental conditions are summarized in Table 2. 
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. Figure 3 


A dry ice pellet suspended in cloudy air, in the FROST tunnel. 


is clearly visible 


The wake 
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Temperature Velocity LWC Time To Initial 
Run No. (°C) (in H20) (ms*') ( gm"9) Breakoff Diameter Comments 
Before After Ap (sec ) (cm) 
1 26.5 26.5 0.10 7 0 40 1.27 
2 26.5 26.5 0.40 14 0) 70 1.27. photos every 10 s 
3 2665972625.> 12707 28 0 100 1.27 
L 18.0. 18.0, 0.43, 14 0 228 1.48 
5 18.0 19.5 0.45 14 0 300 1.54 photos every 10s 
6 19.5 19.5 0.43 14 omen 290 1.53 
7 1955 99:5 0.43 th = 170 1.52 
8 -10.0 -10.0 0.40 13 0 362 1.54 
g -10.0 -10.0 0.40 13 0 396 1.59 photos every 15 s 
10 -10.0 -10.0 0.40 13 0.45 391 1.53 
11 -10.0 -10.0 0.40 13 0.38 300 1.61 
72 -10.0 -10.0 0.40 13 1.40 342 1.54 
13 -10.0 -10.0 0.40 13 1.72 363 1.64 
14 -10.0 -10.0 0.40 13 0 308 1.43 
15 -12.0 -12.0 0.45 13 0 330 1.40 
16 -12.0 -12.0 0.45 13 0 390 1.41 photos every 15s 
17 -11.0 -11.0 0.40 13 0.42 420 1.36 
18 -11.0 -11.0 0.40 13 0.65 402 1-35 
19 24.0 24.0 0.41 14 0 343 Ve52 
20 32651 32605 Oe? als 0 238 1.40 photos every 15 s 
iI 325532. 50. 42a0 14 0 295 1.44 
22 32.5 32.5 0.42 14 3.08 204 1.62 
23 32.5 32.5 0.42 14 3.08 154 1-61 
24 34.0 34.0 0.43 14 0.46 223 1.60 
25 34.0 34.0 0.43 14 0.46 159 1.65 
26 27-0 30.0 0.43 14 0 243 1.32 
27 30.0 32.0 0.52 16 0 247 1.39 
28 32.0 34.1 0.52 16 0 182 1.24 photos every 15 s 
29 34.1 35.5 0.52 16 0 PU? 1.62 for exp 26-30 
30 365,536.06 O.62 0016 0 240 1.54 
31 37503725" 2200" 31 0 80 1.45 photos every 7.5s 
2 sWlos) sWleb rote) sh 0 165 1.52 for exp 31-32 
Table 2 


Summary of Experimental Conditions 
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20 
3.3 Experimental Observations at High Temperatures 


A total of twenty one experiments were carried out at temperatures 
above 0 °C. Six of these experiments were conducted with the sprays on, 
in order to investigate the effect of cloud on the sublimation rate of 
dry ice. Refrigeration problems experienced after conducting 
experiments 6 and 7 prevented the measurement of liquid water contents. 
A rotating clamp was designed for experiments 26 through 32 to eliminate 


the wedge shapes of the earlier experiments. 


The results of the experiments carried out in high temperatures 
with the sprays off are presented in Figures 4 through 12. The smooth 
curves represent the model predictions of the change in equivalent 
circular diameter with time. The jagged lines join the exper imental 
observations obtained from the analysis of the photographs’ recorded 
every 7.5, 10 or 15 seconds. The specific time intervals are indicated 
in Table 2. The temperatures in the working section of the tunnel 
Varied ofrom. “18.0 °C stow 37.55.°C. Velocities in the tunnel were 
maintained close to 14 ms~! with some running at 16ms~'. Experiments 31 
and 32 were conducted at a much higher velocity of 31 ms~! in order to 


test the response of the model to such extreme conditions. 


The differences in the initial pellet diameter, between experiments 
conducted under the same conditions, sometimes varied by as much as 0.11 
cm. Hence, after the photographs were analyzed, the experimental curves 
were shifted to a common initial diameter, for comparison with the model 
predictions. This technique also provided additional information on the 


measurement variability encountered among these experiments. 


eoutarsyed is to belies Sy Biron vagus wy 


;fO ave ide ont Adiv Petculmigg: abu pees 86 


“efesitep cw vi 1} Fe- 5 insite a 5G 5 ters neviari@e KM 


eleven] ase? "se f depen been eae td fe eet uy or seats boat arnien ¥ 


eftremitedss gated ,wfses!h sal lsc feisini ae «] esonsadtih eat 


it. on deem 26 yd belay canihdieee .anotd lone gmne ont—“Tebrat had arbres 


teow ets itiw nati vsqumo 1) 4 Ttem| ib lwlainl noms « of test ide samy - 


10 fo nodgemipth! ‘cools obs beelvon vele supindosy aly ..enelsothaay, 


ia 


=a : ~ " 
zamtersqniot pit os inotaen ged ) 


o 


~ gaolteniftie et ie bute > to tose ald 


16375 ssorp! 19qK9 are Ribs 1 


(73 O32 7 tooo wi? Gi zy revnal 1QAE ‘<i> sro gh apt ae 5 
: ; 5 ev be 7 
-Esnel wae or of 3 


é aah 


x2) Sole AL 3b) “RSP Ane aa pe ’ 


ant <5! dp aE) . sid VM DePrae te ype 4 i 

ponents mt? Vo enohtgheete Jawa ota 18 } | 

* wes ho] eet! tepoey off sul? bathe: ——— 

iqevatste wa oe. ete lee One ware pantiardg at 

s%s efaevesial oma oi Lissase Rail ~via) 0808 4 
mir Oo Tose prisrae one, HI ze eg Ey aT male 

teorws sty af eatetewtey © <a" ak a oN Bt ett aalvat 


= 


> nn} '“en 16 Jo gtisoteyv 446%qit dain 2-46 srasbese time” t¢ bee 


Tale 


gnats (ASS shel ias Asie w isten “63 4 senger was tees” 


% 


inteornl veges on! .neryisgve sou edastpevons of2 wits ,sokt vas 


4 : —_ : = - a wee 
~ Cereal eqs sani” Evers oo vafquobns ya) tigeT ey FADS wena 
_ L re) 


— ¥ : ny i 


_ \ cs - é 4 
~ 


Ze 


The model predictions generally were within 10%-15% of the 
experimental observations . Good agreement was shown for the first 90 
to 120 seconds; however, there was a tendency for the model. to 
overpredict the sublimation rate especially at long times. 
Discrepancies may be due to erroneous estimation of diameter due to rime 
observed forming on the pellet, which was subsequently blown off. Also 
the pellets after some time ,would sometimes slip in the clamp and not 
be perpendicular to the airstream, again giving an over-estimate of the 
diameter. The shape change from circular to wedge also contributed to 
the discrepancies. Another possible source of error may be due to the 
clamp's effect on the air flow pattern around the pellet. The 
variability in actual pellet density could also be a source of 
experimental error. Some pellets could be more dense than 1.4 gcm73, 
and would not sublime as quickly as predicted. A possibility also 
exists that the model assumption of a dry ice surface temperature of 
-100 °C (Fukuta, et al, 1971) may be an overestimate. A colder surface 


would imply a bigger heat flux and hence a faster sublimation rate. 


The purpose in using a rotating clamp in the last series of 
exper iments, was to try to eliminate the shape change problem identified 
in the earlier experiments. This modification did indeed make _ the 
analysis of the photographs easier. The variability between exper iments 
is not as great, as seen in Figure 12, in fact the observations were all 
within 10% of each other. This agreement between experiments could also 
be attributable to amore un i form sample of dry ice obtained from the 


manufacturing company. 
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The results of the experiments carried out in warm temperatures 
with the sprays on are presented in Figures 13 through 15. Comparisons 
with the model predictions for the first 2 minutes show reasonably good 
agreement. The pneumatic atomizing nozzles were not functioning 
properly during experiment 6, resulting in a poor distribution of the 
spray in the working section. This may account for the poor agreement 
with the model predictions. Again a layer of frost was observed forming 
on the surface of the pellet, possibly accounting for discrepancies in 
diameters. The results in Figures 14 and 15 suggest that the liquid 


water content has little effect on the change in diameter. 


A direct comparison between the model predictions and exper imental 
observations when the sprays are on and when they are off is difficult 
to make because the experiments did not begin at the same diameter. 
However, by shifting the curves to a common diameter a comparison can be 
made. A decrease of 38% in diameter iis predicted after 90 seconds, 
(Figure 6), by the model when the sprays are off. A decrease of 64% in 
diameter is predicted after 90 seconds, (Figure 14), by the model when 
the sprays are on. Thus the effect of the cloud is to enhance the heat 


transfer and thereby to accelerate the rate of sublimation. 
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DIAMETER (CM.) 


30 60 90 120 150. ~—«+180 210 
TIME (SECONDS) 


Figure 5 


Same as Figure 4. Conditions for experiment 19 were; 14 ms7!, 24 °C. 
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Figure 6 


Same as Figure 4. Conditions for experiments 20 and 21 were; 14 ms-', 
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DIAMETER (CM.) 


30 60 90 120 150 180 210 
TIME (SECONDS) 


Figure 7 


Same as Figure 4. Conditions for experiment 26 were; 14 ms~!, 28.5 °C. 
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Figure 8 


Same as Figure 4. Conditions for experiment 27 were; 16 ms~! ,. 31.0 °C. 
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Figure 9 


Same as Figure 4. Conditions for experiment 28 were; 16 ms~', 33.0 °C. 
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Figure 10 


Same as Figure 4. Conditions for experiment 29 were; 16 ms- |". 34nb 1G 
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Figure 11 


28 


Same as Figure 4. Conditions for experiment 30 were; 16 ms~!, 35.5 °C. 
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Figure 12 


Same as Figure 4. Conditions for experiments 31 and 32 were; 31 ms-1, 


S7o5 °C. 
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DIAMETER (CM.) 


oA 30 90 120 150 180 


TIME (SECONDS) 
Figure 13 


A camparison of theory and experiment at warm _temper atures with the 
sprays on. The heavy lines’ are the theoretical prediction, and the 
jagged lines are the experimental results. The conditions for 
experiments 6 and 7 were: 14 ms~!, 19.5 °C. LWC was not available. 
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Figure 14 


Same as Figure 13. Conditions for experiments 22 and 23 were; 14 ms7!, 
32.5 °C, 3.08 gm-3. 
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Figure 15 


Same as Figure 13. Conditions for experiments 24 and 25 were; 14 ms-!, 
34.0 °C, 0.46 gm-S. 
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3-4 Experimental Observations at Low Temperatures 


A total of eleven experiments were carried out at temperatures 
below 0 °C. Six of these experiments were conducted with the sprays on, 
in order to investigate the effect of cold clouds on the sublimation 


rate of dry ice. In all cases photographs were taken every 15 seconds. 


The results of experiments conducted at cold temperatures with the 
Sprays off are presented in Figures 16 and 17. All these exper iments 
were carried out at an airspeed of 13 ms~!, and temperatures of -10 °C 
or -12 °C. The sublimation of the pellets is expected to take longer, 
since the temperature gradient between the environment and the pellet 
surface is not as steep as when the environmental temperature is above 0 
oC: The model predictions show very good agreement with the 
experimental observations. The tendency for the model to overpredict 
the sublimation rate is not as obvious as in the high temperature cases, 


but it is still evident. 


The effect of ambient temperature on the sublimation rate can best 
be seen by comparing Figures 6 and 16. Using the shifting curve 
technique, the model predicts a 10% change in diameter, (Figure 16), 
after 60 seconds of elapsed time, when the temperature is -10 °C and the 
airspeed is 13 ms~!. When the temperature is 32.5 °C and the airspeed 
is similar, the model predicts a 21% change in diameter, (Figure 6), 


after 60 seconds of run time. 


The results of experiments carried out at cold temperatures with 
the sprays on are presented in Figures 18 and 19. Exper iments 
10,11,12,and 13 were conducted at an airspeed of 13 eke Ee temper ature 


of -10 °C, and liquid water contents of 0.45 gm-3 , 0.38 gm~3, 1.40 gm73, 


t aa mt laa | = 


7 > a 


aruiseaaret #6) *6 ered 


2 A . 4 


sesutasenre: ss so. belawes Sse <i Seanad 


fo evoiee wit athe Gee aenep ste 2repet pene sauna Yo xl 


: 
nelseriidee sti no thule tia Fe 178 or aaa ae 


so5k 2! vie aves id ange rn? ty és9a0 i1@ mal wat oe 7 


~ ds = 


itiw entuse see Glos fa balay 2 kom *Sg42 1 a2 tose ot 
j / 
sna eee seeds Tan Vf ker ot es wg hh pP shaban 24 te 


2° Of} 1 aewin ieqred teo , “em £1 Go Deeeevie oe +s sug 
: 

rAgne sins oF bagoeqes 2! at’gt ion SA2 16 pe ow oe c, 

~ ar; wis pee ineano*|vire wt) neoedet srobanre aw XS veges oto & 

8 ¢i s juste segues: tz euppetios | view wif caw Rea GE) oo frou att 

otf wy Ty { alee 2104 bias ‘vw ots a ‘soe ao 


jaiteigteve er fester ate (tot peabter eT Daher fatnen! * 


\@eaer swosvsmet agit att dt ee evueteuoler- teh etaes anideiet Tei 


; -iretiws iFige he PT } 


Nee “ 


ot o os 
Sena feo ote raifag ville at 6s om Je ronst Jub ie 1 soothe oe 7 
= 


J a 


eee paleo he ats aiel ell at tree 2! 2 wet? “pnts: ed varies 


lat « wolts ene bb “it wASim POIs - | bd “oy ‘Iekes ms oto 
. 7 . : sg re : 
ons tre 2° Of= gl Btute eS 1) etki emia bean Yo sess 02) ote 

. Songeie gp baw de eos ci asd caer ee nat : ae bt ah wk ~ 


haere ¢ Slane | & u venue $B  ERottess tabon 7s “pratt 
yeni my Te shee es 


33 


and 1.72gm-3 respectively. The results of the model are in good 
agreement with the observations. As the liquid water content increases, 
the observations show a decrease in the diameter, but they are all 
within 10% of the model prediction, suggesting that the liquid water has 
little effect on the sublimation rate. A decrease of 14.2% in the 
diameter after 90 seconds, (Figure 16) is predicted by the model when 
the sprays are off. A decrease of 15.7%, (Figure 18), is predicted by 
the model when the sprays are on. The inclusion in the model of the 
deposition of water vapor and its latent heat, when the sprays are on, 
appears to have little effect on the sublimation rate at low 


temperatures. 


The behaviour of the numerical model under varying conditions is 
summarized in Figure 20. The ambient temperature has the largest effect 
on the sublimation rate. The warmer the environmental conditions are, 
the more rapid the sublimation rate. The effect of a cloudy environment 
tends to enhance the sublimation rate of dry ice, especially in 
conditions above 0 °C. The conduction—-convection term is the dominant 
heat transfer process in all cases. The latent heat due to sublimative 
heating of water vapor is a supporting tare which is more important to 
take into account in warm enviromments. Its effect in a cold and cloudy 


environment is almost negligible. 


While writing this thesis, it came to my attention that a paper on 
experimental studies of dry ice nucleation (Horn et al.,1982) was 
forthcoming. A subsequent telephone conversation revealed that they had 
questioned the Fukuta surface temperature and used a surface temper ature 
of -78 °C. This assumption was noted earlier in my report as a possible 


source of error. | therefore, decided to rerun some of the numerical 


ee 


bag nl ts eRe _sas- ve atluwe aft 
“yeegeveel -snethos “sai stupit (eet eA eqolgs 


ite oie vans ave , vedsme the Sob ni recy 
| oan: aa%ee Ofori ead Sueg thik sail eatta eee ‘teben, 8 
ais ni Bat Feo sense A (ste soiseed Tee ernie 
netw tshew ets vd neo bed ek tat ay . She sea he pane 
| ra ; ioiteng -@t CEL moped) sFA OT wali a «Kio ove mya 
ex Pe lithe gi mi noiaulow. oT ae aa axenne! i bate iso. 
5 aN “wa gig mete , tas poneel 2) Sra “1cRAy. sb ae tw nelatag 2a 


~~ 


at m 6ate> «6(nebteetaee one ne jastts, #124sil ome! La Bi =™ 


et areifiseco balytev som fetom fas! ange ; wit bh wot vanes alee "a 
goslia f29c7el sf) een siseie tage? fislivs sat 4OF swell fi bai 5. y2 Z 


2% e|fol sy ieaoe Lasegmeiensd me sty \oetved ett ost natzonitt ae od # 
themmowlvns yunlS 6 Io saute st | ore wo} tit ( ci 47 big’ von 
nt vileipedes (46+. von. Vo. ates phil ied: ~ orate oF 2 
jheniagd Sif ~ Fae ed nd (Sogwtoseriols tydnot wf > | 2, pooh. notsinany, 


os 


evicemi (vue of seb Seer Janta! @¢T _.dpeen Pe al dzecowy SETA foot 
: i 


ot Snsrvoqm): syom ei viaidw wee polotegode & £1 Waay ~ St uM Yo oniseed 
x 


vmole bee eo 2 BP tuetts eff .ateiag)vte mien 4! Sapodee G2el pte? 
| 


seiciglican secmte mt Jreminriey f we 


no Fee Ss tant niianesi« ym a2 Sita 3) getesiy 213 caisiw Sti | 
hor NST, le te net) netteniaen sot “Wit Ne hl ings - hiSeweth vem 
min yore = peteow a> t2ntihagelon ia inepanda A «Ort bmats wo? ea 


= sts rag EI a been bas setae Brokat ote a 
scien ae xen ee om nobel 4 


34 
exper iments with a surface temperature of —78 °C. 


Results from the 4 different conditions (warm and dry; warm and 
Satur ated; cold and dry; cold and saturated) ane presented in Figures 
21 to 24. The model predicts a slower sublimation rate for all 
conditions considered. These predictions show better agreement with the 
exper imental observations, confirming the earlier concerns about’ the 


surface temperature. 
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Figure 22 


Same as Figure 15. Theoretical prediction determined using a dry ice 


surface temperature of -—78°C. 


42 


aunjzesoduay aoejuns 294 Aap e Buysn 


ove Olde Ost Ost 


poujwsayzep uo} yd|pasd (ed) Ja4054L 


€z aanbl4 


(SGNO53S) 3Wwil 
ord 06 09 0 


€ 


*I08/- 40 
*Q, a4nby4 se awes 


) 8°O 


(WD) 4aLawvid 


fF steph 


‘eudere@es sestwe got (vt 4 gilau Hanimiessh nolsalhe yy [sol seeegtt «0! ayvet?. ce ome, 


PR- Ye" 


t 


sé 


43 


aunqessduay 


soe juns 904 Asp e Buysn poujwacqep uo ol paad [edi yasoayy 


#Z e4nbI4 


(SGNO93S) IWIL 
ost OS} Ozt 06 09  O€ 


°2,82- 40 
QL e4nbyy se owes 


(‘WD) YaL3WVvIG 


eE opty 


gouge sane: scb Te aol yh 6 wad wa hast lands! nat do ibaty inoitssoedl Bt aah ie >a 
| Et= Yo “ 
. v — 


CHAPTER 4 


THE FREE-FALL EXPERIMENTS 


4.1 Free-Fall Model Description 


The theory developed in Chapter 2, was incorporated into a 
free-fall model. In comparing the predictions of equation 10 with the 
wind tunnel experiments, the wind speed which enters Re and Nu is taken 
to be a constant. However, for calculating sublimation during 


free-fall, the wind speed must be the pellet's terminal velocity. 


An object or body in free fall attains its terminal velocity when 
the drag and buoyancy forces are balanced by the gravitational force. 
The development of the expression for the terminal velocity of a 
cylindrical object falling freely in air, in a horizontal orientation, 


is given in Appendix E. The terminal velocity is given by: 


1/2 


‘iO! GD 
ve =f =) (12) 
Toc oe 


where g is the acceleration due to gravity, p, is the air density and 


Cy is the drag coefficient of the cylinder. 


Drag coefficients for objects of given shapes depend primarily on 
the Reynolds number. The drag coefficients for circular cylinders of 
several different aspect ratios are presented in the following table. 


All the values of Cy are valid in the Reynolds range of 10° <Rei< 10°, 


(Hoerner, 1965). 


4k 


wy yigiw Of, Heizeume %6) 2neltoitem, “Sls eal was. A) wba 


ede? ¢1 so ome OF 27eges Asie ly vemie biti aig pot a il 


* RITANQ ~ 


| SPARKS 18 RR BR 


s cont bite reqiecd! ten of mapa “aj! todabeinn: el 


~ 2 


a 


Hie @ he aug igthteiaen “ah Sewer i ta & 


Vitlootev tealimed <'saltae sig Si oon dsage we aty ,! 
SO 
wty yitoulee’ Tentunes 27! entedan Lhe evr nel via: 7 ——r \ 


“hy 


sO5O7 feMol sel i ussp ons yd Spare lA tetw penta Money ous tens wie : 
to. &ateates fenieiee mete . ea ih exengas wet 3 | 
nelses atva fears] ne 3 af 06 A) yee oe Aza. 
vil oor | Be). inv leniwres sat Seeeraeahts 


brie Staten te otP al GO yy fives x Sib nels £6 VST ot 


—< 


y 


ovebnifys s47 Yo tes! at Ttaup 


qa . 


BO gFinnney, brageb gsagts wo tw aii 403 aren OTR 


Yo emt Fre wivotio ww? esnat sien He: sa ware 
oidas priweiigh sis ni weacom 18 noltey sadhen 3 


_ 
Bie “Of Wo ep, ani al ny 35 oe 


45 


OL SS SS -S --<-- l S- sSSSSS 


Aspect Ratio Drag coefficient Reynolds No. range 


L/d 


1 0.63 10°< Re < 10° 
> 0.80 
10 0.83 
20 0.93 
30 1.0 
za 2 
Table 3 


The Drag Coefficient of Truncated Circular Cylinders 


A drag coefficient value of 0.8, corresponding to an aspect ratio 
of 5 was used in the free-fall model. An atmospheric lapse rate of 6.6 
°km7! was used in order to specify the variation of the ambient air 
temperature with height AUninG the fall of the pellet. The change in 


the aspect ratio due to sublimation was ignored. 


4.2 The Free—Fall Experiments 


A series of free-fall experiments was conducted, in which 
cylindrical dry ice pellets were dropped from an aircraft flying at 305 
m, 460 m, and 610 m above an unused runway (915 m MSL). The intention 


was to measure the time taken for the cylindrical pellets to reach the 
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ground, and the final diameters they achieved, and to compare these with 


the free-fall model predictions. 


Prior to the drops, plastic bags with about 1-2 kg of pellets 
(cylinders with average diameters of 1.27 cm, varying between 1.2-1.5 
cm, and average lengths of 5.0 cm, varying between 1.3-6.8 cm), were 
prepared. As the aircraft approached the runway, the contents of a bag 
were dropped through an opening in the bottom of the plane. An observer 
on the ground followed the aircraft through a pair of binoculars, and 


timed the fall of the pellets with a stop watch. 


The model predictions of the time taken for the pellets to reach 
the ground, and the final diameters they would achieve are presented in 
Figure 25. The solid lines show the pellet diameter as a function of 
height, for the three different release heights. Superimposed on these 
lines are isochrones, in seconds, represented by the dashed lines. For 
pellets of initial diameter 1.50 cm, the predicted fall times for the 
three release heights are 15.7, 23.4, and 30.1 seconds, respectively, 
with corresponding final diameters of 1.42, 1.40, and 1.34 cm. The 
observed fall times for the first pellets reaching the ground were 18, 


26, and 36 seconds. 


The difference between the predicted and observed fall times could 
arise if the pellets encounter thermals as they fall, thereby retarding 
their fall. Another reason for the difference could be that the drag 
coefficient on the pellets is somewhat larger than’ the value 0.8 
considered in the model prediction. If a higher value of the drag 
coefficient was used, the model would then predict fall times closer to 


the observed fall times. There appears to be a larger discrepancy 
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between the predicted and observed fall times as the release height is 
increased. Scintillations observed through binoculars as the pellets 
fell suggest that tumbling of some sort was most probably occuring. The 
pellets in this case would achieve a higher velocity than they would 


with the longitudinal axis horizontal. 


An attempt to compare final diameters between model and exper iment 
was not successful. The pellets dispersed considerably in the 
horizontal as they fell, (also possibly indicative of tumbling) and many 
were lost in the tall grass bordering the runway. The next time such 
experiments are attempted, it would be advantageous to find an area with 
grass that has been cut. Those that were found on the runway were very 
much smaller than the predicted sizes. It is suspected that the impact 

at about 17 ms~!, on a hard runway surface probably shattered the 
pellets, rendering invalid any final comparison of diameters. An 
experiment was conducted to determine if such was the case. Some dry 
ice pellets were dropped off the roof of a building (height of 15 Mm). 
The impact at abet qi ms~! on the pavement below, indeed shattered the 


pellets. 
4.3 A Free Fall Simulation 


A series of numerical experiments were conducted by using the model 
described in section 4.1. This simulation was conducted to determine 
the amount of dry ice that sublimes within the supercooled region of a 


cloud. 
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Dry ice pellets with initial diameters varying from 0.4 cm to 1.0 
cm were introduced at the -10°C level. The emulacion also assumed an 
atmosphere with no vertical motions. The results of these numerical 
experiments are presented in Figure 26. The solid lines are the pellet 
diameter as a function of height or temperature. The dashed lines 
represent isochrones in seconds. The results indicate that pellets with 
diameters between 0.5 cm and 0.6 cm are suggested for efficient use of 
dry ice when dropped from the -10 °C level. The term ‘efficient! means 
the optimum size of dry ice pellet that will sublime within the 


supercooled region of a cloud. 


The dry ice pellets would have a longer residence time in the 
supercooled region of a cloud with wupdrafts. In fact, they might 
‘completely sublime above the 0 °C level in a cloud with a vigourous 


updraft, indicating that larger diameters would be required. 


ee ae 
Gel og a Oe. mon’ ny everenmi ietaial athe 


ne bemeoe ua ts roisk time’ ot Jn tive! Pot~ wa eee 


ce 

tant semen seata to edtuaed WAT .enctien bap thahiy ne ee 
| a eee 

sol(feq otf 21s cent! of foe Sm. pe a mught a ‘neo iacng Le | hs 


oe 


eoal!t bedzad rT 2 36 rHqme? We, 24g. ro nmiane 80% 


Agtw avfelleq seas adel ml stfuzes af . &batecrae a ssnavddae| 4 


us: Sauloltts 16) bePeeton? 2 4.9 Of & 2. ve fain: 
> : f , F 7 = oll 
orges ‘tneioltie! wae eff <Tewet 3° 0 be ont vet pedehont Reldte aat ; 


ait aleghe smttde (lle sede SPT itp. GSP yo ‘a oxte ‘aml 


"> Say 


wots 6 to aelpe te 
Fi . oh | 
oi? Al pe " eorah fee ‘aga! = B sve pluie | emot eet ir. Mt 


y 
Solm. vety . foe? in 27achy «6 Senha é ¥9 , ele te 


a , 


woweonts & dilv Quelp « «i fevel a* 0 ete aved oral due hase tque 
haar take 3d tidus etapa |b. mele! 7aRns enti 1 ei 


49 


*Spuooas 
ul} SaUuO‘YyDOS| SY dae Sout, peysep sy *sqyBbyay aseayos aa4yQ 404 BYHyay JO uoyzOuNJ e& Se 
yajoue ip yo_ped ayy a4e sous, PILOS SyL *squaw}sadxa | [ejJ—-994J BYR JOJ SUOJZOIpasd | apow 

Gz eanb1 4 


(WD) waLlawvid 
s*t vt et 


Oot 


o 
3 
(W) GNNOYD JAO#V LHOISH 


° 
8 


009 


| . et wtuel 4 : e 
FFs —e feliSq séx ew-eudi! bifon. off tntodigges . tie t-2n17. ods oer ast Feyticots G. = 
: 6) eandsdood! sz. ors aanil bedeels eat .etdeied peagher. sows ~ = Yo solzonut © 2s ) 


| viacsiee iin Me | . <abaiocain “eS 


= { 


50 


*spuooas ul Ssauosyoos| ade sauit pe 
jo uoyzOuns e se sa owe |P qajpad ayz ase sauyt PLLOS UL 


gz 34nb} 4 


90 sO 


o°t 6°0 8°0 L°0 


(D:) SuNivdsadwal 


ysep ayL sounqesoduay 40 3ybyay 
*Bulpeas adj Asp Jo uoy we_Nwys V 


cost 
o0zz 
1 0092 
000 
OOvE 


oost 


(W) LHOISH 


002 


oo9v 
ooos 


OcvS 


ooss 


of saugi4 . é* ‘ 
“*e Neal lati? ¢ ee re Fane sh talieg =tt aie 2401? Bites sfT pi imese eal Vb To oO 
 -8bPooee OF tenstfaeet ots cent! heiteeb ant »e nem 


CHAPTER 5 


SUMMARY CONCLUSIONS AND RECOMMENDATIONS 


The objective of this thesis was to utilize both heat transfer 
theory and experimental observation in order to examine the sublimation 


rate of cylindrical dry ice pellets in clear air and simulated cloud. 


A numerical model assuming a heat flux balance at the surface of a 
cylindrical dry ice pellet was developed. The sublimation was assumed 
to occur from the cylindrical surface only and not from the ends of the 
pellet. It was also assumed that the sublimation took place uniformly 
both around the circumference and along the length of the cylinder, so 
as to retain its cylindrical shape. The airstream was assumed saturated 
when simulating a cloud, and was perfectly dry at all other times. The 


contribution of cloud droplets to the heat transfer was ignored. 


A series of thirty two laboratory experiments were conducted in the 
University of Alberta's FROST icing-wind tunnel facility, to test the 
theoretical model. These experiments were carried out in both cold and 
warm environments. Some of the experiments were carried out with the 
sprays on in order to determine whether a simulated 'cloudy' environment 


had any effect on the sublimation rate. 


The ambient temperature was found to have the largest effect on the 
sublimation rate of the dry ice pellets. The warmer the environmental 
conditions, the more rapid was the sublimation rate. The effect of a 
saturated environment also enhanced the sublimation rate. This effect 
was more pronounced at warm-temperatures, and almost negligible at cold 
temper atures. In spite of the use of several simplifying assumptions, 


the cylindrical model, (using a surface temperature of -78°C), predicted 
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the sublimation rates of dry ice pellets to within 10%, when compared 
with wind tunnel observations. The model did not perform as well using 


Fukuta's estimate of the surface temperature of -—100 °C. 


The theoretical model was then modified to numerically simulate the 
sublimation rate of cylindrical dry ice pellets in free-fall. A series 
of nine field experiments was carried out whereby dry i was released 
from an aircraft flying at specified altitudes, to verify the model 


predictions of fall times and final diameters at the ground. 


The model predictions of the fall times were within 10-15% of the 
observed times. Attempts ‘to compare final diameters with the model 
predictions were not successful. The impact on the runway surface 
shattered the pellets rendering invalid any final conparisons of 


diameter. 


A numerical simulation of a dry ice seeding experiment was 
undertaken, to determine the amount of dry ice that sublimes within the 
supercooled region of a cloud. The results indicated that pellets with 
diameters between 0.5 cm and 0.6 cm were suggested for efficient use of 


dry ice when dropped from the -10 °C level. 


The FROST icing wind tunnel was an appropriate facility for 
conducting the sublimation rate experiments. This facility may also be 


useful for investigations into ice crystal production rates. 
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APPENDIX A 
Determination of the Nusselt Numbers from 


Experimental Observations 


A method of obtaining a relationship between the Nusselt and 
Reynolds numbers directly from experimental observations was 
explored. Equation 10 applied to the dry conditions’ was 
re-arranged as: 

eel aan at (1) 
as eae 

A series of four experiments conducted in dry air (8, 9, 14, 
and 19) was used to calculate the Nusselt number. A three point 
averaging filter was applied to the observations of dD/dt ys) 
smooth out the large variations. The resultant smoothed values 
were then inserted in the above expression and a Nusselt number 
calculated. Corresponding Reynolds numbers computed from the 


experimental data, were substituted into the Zukauskas expression: 


Nu = 0.26 Re°’® (2) 


The results of the calculations are shown in Figure 27 in’ the 
form of a scatter diagram. The Nusselt number determined using the 
Zukauskas expression is plotted on the abscissa, while the Nusselt 


number obtained from experiments is plotted on the ordinate. 
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The 1:1 line is superimposed on the data. A linear regression 
analysis was performed on the data using the method of least 


squares. A correlation coefficient of 0.7 was calculated. 


The Nusselt calculated from the experiments seems smaller than 
the Nusselt from the dZukauskas. This is consistent with the 
observation that the sublimation rate experimentally is generally 
lower than that predicted by the model. Errors in estimating the 
change in diameters explain why the scatter seems largest at _ the 
largest Nu. These results suggest that the observational data are 


not good enough to derive Nu. 
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Figure 27 


Scatter Diagram of Nusselt Calculations 
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APPENDIX B 


Derivation of the Q. Term 
es 


The total latent heat transfer due to sublimation of water 


vapor onto the dry ice can be expressed as: 


lee ete pee) (1) 


where h is the mass transfer coefficient, and P. and 2, 


are the 
water vapor densities (absolute humidity) in the ambient 


environment and at the surface of the pellet. 
The equation of state for dry air can be expressed as: 


P, = oR T 2) 
where R' is the specific gas constant for dry air. 


Water vapor in the atmosphere behaves as an ideal gas to a 


good approximation and its equation of state can be written as: 
rie. (3) 
where e is the partial pressure of the water vapor, Pp, is the vapor 


density and Ris the specific gas constant for water vapor. 


The total air pressure P, according to Dalton's Law 


is equal 
to the 


sum of the partial pressures of the dry air and the water 


vapor. 
peace (4) 
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Then the equation of state for dry air can be re-written as: 


= (5) 


Considering the ratio of Q/p: 


le) 


ee = e= (6) 
0 


or 


Pp =p E= (7) 


The mass transfer coefficient he in equation 1 can be given 


heh 6 
D 


Wa 


Sh = 


(8) 


where Sh is the Sherwood number, and Dis the diffusivity of water 


vapor in air. 


Using the result that the vapor density can be related to 
vapor pressure according to equation 7, equation 1 can be 


re-written ass: 


oe kl, mOL(e.-e.) (9) 
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Substituting for Bis according to equation 8, the mass transfer 


coefficient kK can be expressed as: 
= -1 EO 
k Sh ee D 5 : (10) 
or 
Ke Dee Dae ashe Rua ato. (11) 
: R! 
since € = aa 
Vv 
then the mass transfer coefficient may be expressed as: 


k = ae <1 rag 12 
(Oe Y Dy MeeSh REY oF (12) 


Now consider: 


eT S ‘& 
Kee am ones belae 


ep a ts be ch 1 
h k. D7! Nu (13) 
a 
aa! =-l] 
Bn Gaal (ai (14) 
h Kk, aa} 


since Sh = cRe™ Sc" 3 and Nu = cRe™ pr” » where Sc is the 
Schmidt number defined as V/Dwva ; and Pr is the Prandtl number 


defined as v/K . hences 


; Pua Fy 1 (Se) : tas 


LS 
Hae Pr 


however k =KOC « Where « is the thermal diffusivity. 
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Then equation 15 can be re-written as: 
whl 
Sn (i) as es) 
h K (a Pr 
R Pp 
Kn (Pua) (se) " wer 
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(16) 


(17) 


(18) 
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Hence, using a value of 0.37 for the exponent n (Morgan, 


1975), equation 9 can finally be written as: 


p 63 €] 
= cake mnaee - ; 
QF. roun( £2) Pe a) value) 


(19) 
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APPENDIX C 


An Order of Magnitude Analysis of the 


Heat Balance Equation 


A numerical example is presented and the 


relative magnitudes 


of the Qu Q , Q. and Qy terms are evaluated. The following 


conditions have been used 


D 


t 
a 


Re 


Nu 


es r 


1.5x10-°m L = 5.0x10°m 
2026 t.= -100°C 
5.67x10 ° Wm K~ 7 = 3.14159 
k,NuD Wm K7? 

Deus x1Oie ex Ole btn) Wm? K+ 
0.26Re™® 

Vv 
Vem moce 

1.718x10 ° +5.1x107” (t_) kgm” $7 
2.738x10° kgm / si; 
2.147x107> m’s ° 
9.779x10° 

64.4 

3.89x107 Wm? K7? 


1.67x102 Wm7*K7? 


The heat transfer due to conduction 


evaluated as: 


Q = mDLh(t.-t,) 


= -47.2 Watts 
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in the calculations: 


Vv = 14 ms? 


p= 12275 kgm ° 


and convection is 
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The heat transfer due to sublimation of water vapor 


dry ice, is evaluated as: 
63 ] 


Be een 
Q., = -roLn(£) Few (eq (tg)-e, (t4)) 


where € = 0.622 
Pr = 0.711 
Sc = 0.595 
35 = 1.005x10° Jkg? K7 
a 6 -1 
V2 = 2.823x10° Jkg 
p = 930 mb 
e,(t,) = 23.373 mb 


Qog= -19.3 Watts 


The heat transfer due to radiation is evaluated as: 
= - 4 2 4 
0 Oo Tame © Thea eenDE 


where the emissivities are assumed to be 1 


ue = 293.15 K 
Te Re y/ Sy Uy a 
Q = - 0.867 Watts 


onto 


The heat transfer due to the droplets is evaluated as: 
=- c 0 0 Ee 
Q4 (c R(t. Ones, R Le - oi R (0 C P2)) 


where c= 419x102 Jkg7 K-? 


1.74x10° Jkg-!k-2 


oO 
i} 


Le = 3.34x10° Jkgm* 


Rw 


defined as: EVwDL 


is the mass flux of droplets collected per unit time, 


and 
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Where E is the collection efficiency, assumed to be 1 (this is 
an over-estimate. A precise value cannot be obtained because the 
droplet size spectrum was not measured.), V is the velocity, and w 
is the liquid water content. A value typical for a small cumulus 


cloud is of the order 1.0x107? kgm™’. 


Therefore R, = 1.05x10° kgs"). 


In the term each component is evaluated as 


a = 0.88 Watts 
b = 3.51 Watts 
GC = 1.83 Watts 


Qy= -6.22 Watts 


It can be seen that radiation only contributes at the most a 
small percentage to the total heat exchange and may be ignored. 
The heat due to the droplets is not a small amount as can be _ seen 
by the calculations. The contribution to the heat balance, though, 
is unknown since we do not know at what stage in the cooling 
process the droplets leave the surface. In addition, during the 
experiments ice accretion was not observed on the surface of the 
pellet. It is for these reasons that the heat contribution due to 


the droplets is ignored. 


WO PEE hs iy i 
: 8 ; i a) 


sT otat) i st-@? bonvase svotahol tre nataeeriin 
- o wi 
849 seunaed wantesed) oo torres satel a2i meng # 


a wn ~“tiseisy Say a? 


eviimee> tiaonmk @ wT tsolave-s 


—_ 


e 


6 tion Sy Jn sesuci vtew> sf os es Tel eid. rege vind - 
wosnrp! ed vein Bap yn oO RES? Peon bebe At os: anti thn 
Mae 6 Hen 4 Pease ae & Jes e! stela et ang Ne “eb eed tT 
, Apuerts iaihade 4 tld og y WU tadINIGS eat | Bis sricvie 
Bal lees awit ‘ab sper fede 16. pert for obs ‘ey eonihe’ nn a 
sity pstcwt eielthebe a! ssdatebs @lPieyant | ety (Gort wie vi missy 


omy %% @oetoup SH fo erent sonieae Apizosrce ae 2 


oo a relat jenn Aeon -a19 Jeg enneesy seis me al wi 


spans. al 


APPENDIX D 
Analytic Solution to 


Sublimation Rate Expression 


The sublimation rate expression (equation 10), is written as: 


RES 
dD _ a [te + Hae) ane ae | (1) 
t De beD Ss sa Sc ae a) ce 


the terms within the square brackets are a constant with time, say 


-C. Then: 


See Ce (2) 


Substituting for the Reynolds and Nusselt numbers, 


_ vd 
Bea io Nu = 0.26 Re™® (3) 


Then: ae. te 
dD 2k V D 0.26 
mrs Or C 
dt oer aoy 
(4) 
dD _ = Taian 
1E 
wheres 
0.6 
2k V 0.26 (5) 
K = c 
ep Ee 6 
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which is also a constant in the wind tunnel, but not in_ the 


free fall experiments. 


This expression may be integrated to give: 
D =n(Dih* fo taiwhke) elas (6) 


where D; is the initial diameter. 
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APPENDIX E 
Derivation of the Terminal Velocity 


for a Cylindrical Object 


The forces acting on a cylindrical object falling freely in 


air are depicted in the following schematic. 


Z Vv 
(+'ve ) 


The equation of motion for this cylindrical object, negiecting 


buoyancy which is expected to be small for dry ice pellets, can be 


written as: 


(1) 


where m is the mass of the cylinder, g is the acceleration of 
gravity, and ce is the drag force. The negative sign appears on 
the left hand side of this expression because V is considered to be 


negative when directed downwards. 
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The drag force (Huschke,1959) can be written as: 
2 
be 2 
F) = 1/2 ACyP.V (2) 


where A, is the cross-sectional area of the body normal to’ the 


airstream. For acylinder this is: 
Aca Det (3) 


where D is the diameter of the cylinder, and L is the length. 


Hence,the drag force acting on a freely falling cylinder, 


oriented with its axis horizontal, can be written as: 


Fy = 1/2 DLC pV" (4) 


The mass of the cylinder is: 


2 
GheoPaate (5) 


Substituting the expressions for the drag force and the mass 


of the cylinder, into the equation of motion, we obtain: 


(6) 
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where 


pa (7) 


The terminal velocity V=V_ occurs when the net force acting on 


the cylinder is zero. 


hence 


: (8) 


Or 


oo (9) 


The terminal velocity for a freely falling cylinder can 


therefore be expressed as: 


’ 4) (10) 
T 2) 
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APPENDIX F 


Computer Program Listings 


This appendix contains the FORTRAN listings for the computer routines 


used in and developed for this thesis. 
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SUBLIME .FOR 


2 eee ieee a I fe tc RIK kf Kk Rok io tk ok kaka aR ok ok ook ko ok oe oe ee 


THE NUMERICAL SIMULATION OF DRY ICE SUBLIMATION USING 
HEAT TRANSFER RELATIONSHIPS. 


THIS MODEL IS VALID FOR CYLINDRICAL PARTICLES. 


We HE Ae eK tee fe ic oe oe oe oe fe ic ei oko ke ke oie ke ok ok io eo ok oko ok ok ak ak ak ok ak ok ie keke ie ai ak kc ake teak ok ae ak ok ok oe ok kok ak ok ke ok ok 


NNgagNNnNNANANANNANAANAANANNO 


REAL NMASS,NRAD,NUSELT,MBCON, INCON,MU,NU,MASS 
REAL NDIAM,LHS,KA,NUMR,KAPPA,LHW 


DEFINE SOME CONSTANTS 


aaa 


RHOCO2=1.40 
CM=2.54 
MBCON=1000. 
INCON=2.49082E03 
LHS=587E07 
LHW=0.0 
PI=3.14159 

RE=O. 

NUSELT=O. 


OR ee ee ee oe oie eke ee kee eK ee eK eK ke i eK oe ok KK Ko ko oe ok KK ok oe ok ok ke ok OK Ok ke ok 


FORMAT STATEMENTS 


se oie fe ee Re oe eo i oo i Ke oe eof oe ic oe oR eo oe eo ke oe oo eo ie ke ok oe ie ok eK io OK KK oe oe ok oo ok ok ok ok ok ok 


FORMAT( 1X, “INPUT---TIME INT.,CYL.DIA,CYL.LEN. ’) 
FORMAT(1X,’INTRODUCE INITIAL CONDITIONS---TA,PA,DELP,TS,VEL,WET’ ) 
2 FORMAT(’1’,42X, ‘THE NUMERICAL SIMULATION OF DRY ICE SUBLIMATION’ ,/ 
+.) 
16 FORMAT(’ ‘’,33X,’DRY CASE---SPRAYS OFF’,//) 
17 FORMAT(’ ’,33X,’WET CASE---SPRAYS ON’,//) 
18 FORMAT(’ ’,33X,’TIME INTERVAL (DELT)=’,F4.2,’ SECONDS’ ,//) 
13 FORMAT(’ ’,9X,’INITIAL CONDITIONS: ’,4X, ‘PELLET TEMPERATURE: ’,9X, 
#F6..1, 1X; “CELSIUS” ./. 33x; “ORY ICE DENSITY: ”, 14x. 
+F5.2, 1X, ’GM/CM**93’ ,/,33X, ‘CYLINDRICAL RADIUS: ’,11X,F5.2,1X.’CM’,/, 
+33X, ‘CYLINDRICAL LENGTH: ’,11X,F5.2,1X,’CM’,/,33X, ‘CYLINDRICAL DIAM 
+ETER:’,9X,F5.2,1X,’CM’,/,33X, ‘AMBIENT TEMPERATURE: ’,10X,F6. 
+2,1X,’CELSIUS’,/,33X, ‘AMBIENT PRESSURE: ’,13X,F5.1,1X,’MB’,/,33X, 
+’WIND TUNNEL VELOCITY: ’,9X,F5.2,1X, ‘M/SEC’ ) 
15 FORMAT(‘1’,15X,’TIME’, 15X, 
+’/DIAMETER’ ,15X, ‘MASS’, 15X, ’REYNOLDS’ ,15X, ‘NUSSELT’ ,/, 
+14X%22 (SEG 17%, CCM) 17x, CGM) 7 ,./,) 
14 FORMAT(13X,F7.2,15X,F8.6,14X,F6.2,14X,F8.2,14X,F8.2) 
666 FORMAT(2F 10.4) 
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FORMAT(1X,4F20.6) 


INPUT CYLINDRICAL LENGTH AND RADIUS (CM) AND TIME INTERVAL (SEC) 


WRITE(6, 1) 

CALL FREAD(S..’SR:” {DELT, CYLD. CYLL) 
LF (DELT ~EQ: 0.) 2 GO_10 90 

TIME=0.0 

J=0 

LINE=0 


CALCULATE CYLINDRICAL MASS .(CMAS.) 


CRM=CYLD/2. 
NDIAM=CYLD 
CMAS=PI1*( (NDIAM*NDIAM) /4) *CYLL *RHOCO2 


INTRODUCE INITIAL CONDITIONS: TA- AMBIENT TEMP. IN CELSIUS. 
PA- AMBIENT PRESSURE IN MB. 
DELP- DYNAMIC PRESSURE INCHES 
TS- PELLET TEMPERATURE IN CELSIUS 
VEL- VELOCITY IN M/S DETERMINED 
WET- WET/DRY CASE? 


WRITE(6,2) 
CALL: FREAD (CS. (GR>7.,7A, PA, DELP. TS, VEL. WET) 


TBAR=(TA+TS)/2. 
TBK=TBAR+273.16 
TAK=TA+273. 16 
TSK=TS+273.16 
PAT=PA*(0.9862E-03) 


RD=0.287E0O7 
P=PA*MBCON 
RHOA=P/(RD*TAK) 


IF(DELP.EQ.0.) GO TO 100 
V=(2*DELP*(INCON/RHOA) )**(1./2.) 
VM=V/100. 

GO TO 11 

V=VEL*100. 

VM=VEL 


WRITE(7, 12) 

IF(WET.EQ.0.) WRITE(7, 16) 
IF(WET.EQ.1.) WRITE(7,17) 
WRITE(7,18) DELT 


WRITE(7,13) TS,RHOCO2,CRM,CYLL,NDIAM,TA,PA,VM 


WRITE(7, 15) 
LINE=LINE+3 
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WRITE(7,14) TIME,NDIAM,CMAS,RE,NUSELT 
LINE=LINE+1 


TIME=TIME+DELT 
J=TIME 


COMPUTE THE THERMAL CONDUCTIVITY OF AIR. 


KA=2.43E-02+(7.3E-04*TA) 
KA=KA*1.0OEO5 


CP=1.O005E07 : 
CALCULATE KINEMATIC VISCOSITY. 


MU=1.718E-05+(5. 1E-O7*TA) 
MU=MU* 10. 


NU=MU/RHOA 


PR=0.711 
SC=.595 


SEPARATE HEAT TRANSFER TERMS CALCULATED 


QC=(TS-TA)/LHS 


TF(WET.EQ. 0.) GO TO SO 


A POLYNOMIAL EXPRESSION (LOWE,1976) USED 
TO COMPUTE SAT. VAP. PRESS. 


AO=6. 107799961 

Ai=4.436518521E-01 
A2=1.428945805E-02 
A3=2.650648471E-04 
A4=3.031240396E-06 


T=TA 
EATA=AO+T *(A14+T*(A2+T*(A3+A4*T) )) 


S0=676.967 
$1=-0.0680427 
$2=-0.120855E-02 
$3=-0.518213E-05 
$4=-0.249304E-07 


TM=TS 
LHW=SO+TM*(S1+TM*(S2+TM*(S3+S4*TM) ) ) 
WRITE(6,666) EATA,LHW 
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LHW=(LHW/O.2389) *1.0E07 
EPS=0.622 

QE 1=(PR/SC)**(0.63) 
QE2=-(EPS)/(PA*CP) *(LHW/LHS) 
QE3=EATA 

QE=QE 1*QE2*QE3 

GO TO 20 

QE=0.0 


CALCULATE REYNOLDS NUMBER FOR NEW DIAMETER. 


RE=V*NDIAM/NU 


USING THE SUBLIMATION EXPRESSION DERIVED COMPUTE CHANGE 
IN MASS,HENCE CHANGE IN DIAMETER, AFTER TIME DELTA. 


NUSSELT RELATION ACCORDING TO ZUKAUSKAS 
APPROPRIATE FOR REYNOLDS RANGE 10**3-10**5 


NUSELT=0.23*(RE**0O.6) 


CONST=(2*KA*NUSELT ) /RHOCO2 

DELD=( (CONST/NDIAM) *(QC+QE) ) *DELT 

ND IAM=NDIAM+DELD 
NMASS=PI*RHOCO2*CYLL*( (NDIAM*NOIAM)/4) 


IF(NDIAM.LE.O.) GO TO 25 
GO TO 30 
NDIAM=0.0 
NMASS=0.0 


IF((TIME-J).EQ.0) GO TO 200 
GOTO: 24 
WRITE(7,14) TIME,NDIAM,NMASS,RE,NUSELT 


WRITE(8,888) TIME,DELD,QC,QE 


LINE=LINE+1 
IF(LINE.EQ.56) GO TO 31 
IF(NDIAM.EQ.O.) GO TO 10 
TIME=TIME+DELT 

J=TIME 

GO TO 20 

LINE=0 

WRITE(7, 15) 

WRITE(7,14) TIME,NDIAM,NMASS,RE,NUSELT 
LINE=LINE+4 

TIME=TIME+DELT 

J=TIME 

GO TO 20 


STOP 
END 
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CYLDRP.FOR 


THE NUMERICAL SIMULATION OF DRY ICE SUBLIMATION USING 
HEAT TRANSFER RELATIONSHIPS. 


THIS MODEL IS VALID FOR CYLINDRICAL PARTICLES. 
FALLING THRQUGH AN ATMOSPHERE. 


NAANDAANANANANAANANAANAANAAAA 


REAL NMASS,NRAD,NUSELT,MBCON, INCON,MU,NU,MASS 
REAL NDIAM,LHS,KA,NUMR,KAPPA,H,LAPSE,KELVIN 


DEFINE SOME CONSTANTS 


qaqa 


RHOCO2=1.40 
CP=1.O005E07 
TS=-78.0 
CM=2.54 

MBCON= 1000. 
PI=3.14159 
LAPSE=1./152. 
RD=0.278E07 
GRAV=980.616 
CLAPSE=LAPSE/100. 
DENM=CLAPSE*RD 
EXP=GRAV/DENM 
LHS=587E0O7 
KELVIN=273.16 
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Cc 
Cc FORMAT STATEMENTS 


Cc 


(SE Re Mee aR RRA RS Ae AR (RAR A RC RE GRC Oe ee a cee a hoe CO 


Cc 


C 

( : 

1 FORMAT( 1X, ‘INPUT---TIME INT.,CYL.DIA,CYL.LEN. ,DRAG’ ) 

2 FORMAT(1X,’INTRODUCE INITIAL CONDITIONS---TA,PA,H, WET’ ) 

12 FORMAT(’1’,42X,’THE NUMERICAL SIMULATION OF DRY ICE SUBLIMATION’ ,/ 
+5) 

UU FORMATO? 9". 33X97 DRAG COEFFICIENT: ’ FS.2,/) 

13 FORMAT(’ ’,9X,’INITIAL CONDITIONS:’,/,33X,’DRY ICE DENSITY:’, 14x, 


+F51291X 27 GM/CM# #377 /533X, “CYLINDRICAL RADIUS: ’,11X,FS.2,1X,’°CM’./, 
+33X,’CYLINDRICAL LENGTH:’,11X,F5.2,1X,‘CM’,/,33X, ‘CYLINDRICAL DIAM 
+ETER:’,9X,F5.2,1X, ‘CM’ ,/,33X, ‘AMBIENT TEMPERATURE: ’,10X,F6. 

+2. 4X 9 CEUSIUS’ ¥/,33X,” AMBIENT PRESSURE: ’,13X,F5.1,51xX, MB”, /,33X, 

+ INI TUAL TOROPSHEIGHT: “..10X,.FS8.4,1x, - CET MSE)” .///,40X. “TIME”, 
+5X,’TEMP’ ,5X, ‘PRESSURE’ ,5X, ‘HEIGHT’ ,5X, ‘DIAMETER’ ,/,40X, ‘(SEC)’ 
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FORMAT (38X,F7.3,4X,F5.1,6X,F5.1,5X,F8.1,6X,F8.6) 


INPUT CYLINDRICAL LENGTH AND RADIUS (CM) AND TIME INTERVAL (SEC) 


WRITE(6, 1) 
CALL FREAD(5,’4R:’,DELT,CYLD,CYLL , DRAG) 


TF? COELTEQ.0.) GO T0790 
TIME=0.0 


CRM=CYLD/2. 

NDIAM=CYLD 

CLM=CYLL 

CMAS=PI1*( (NDIAM*NDIAM) /4) *CYLL*RHOCO2 


INTRODUCE INITIAL CONDITIONS: TA- AMBIENT TEMP. IN CELSIUS. 
PA- AMBIENT PRESSURE IN MB. 
H- DROP HEIGHT (FT. MSL) 


WRITE(6,2) 
CALL FREAD(5,’4R:’,TA,PA,H,WET) 


TAK=TA+KELVIN 
PAT=PA*(0.9862E-03) 


P=PA*MBCON 
RHOA=P/(RD*TAK) 
WRITE(7,12) 
WRITE(7,777) DRAG 


WRITE(7,13) RHOCO2,CRM,CLM,NDIAM,TA,PA,H 
WRITE(7,14) TIME,TA,PA,H,NDIAM 


COMPUTATION OF TERMINAL VELOCITY,MEAN TEMP. ,MEAN PRESSURE, 
AND DISTANCE FALLEN DURING THE TIME INTERVAL SPECIFIED. 


X=PI*RHOCO2*GRAV*CRM 
Y¥=DRAG*RHOA 
TV=SORT(X/Y) 
DIST=TV*DELT 
H=H-(DIST*O.03281) 
DISTM=DIST/100. 
TNEW=TA+(LAPSE*DISTM) 
TNAB=TNEW+KELVIN 
PNEW=PA*((TNAB/TAK) ) **EXP 
PMEAN=(PA+PNEW)/2. 
TMEAN=(TA+TNEW)/2. 
TAK=TNAB 
TMK=TMEAN+KELVIN 
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120 TA=TNEW 


124 PA=PNEW 
422 P=PA*MBCON 

123 PAT=PA*(0.9862E-03) 

124 TIME=TIME+DELT 

125 RHOA=P/(RD*TAK) 

126 PM=PMEAN*MBCON 

127 RHOAM=PM/ (RD*TMK) 

128 (8 

129 c 

130 c 

131 C----- COMPUTE THE THERMAL CONDUCTIVITY OF AIR. 
132 c 

133 @ 

134 KA=2.43E-02+(7.3E-04*TMEAN) 
135 KA=KA*1.OE05 

136 (3 

137 C----- CALCULATE KINEMATIC VISCOSITY. 
138 Cc 

139 c 

140 MU=1.718E-05+(5.1E£-O7*TMEAN) 
1414 MU=MU* 10. 

142 (3 

143 . NU=MU/RHOAM 

144 Cc 

145 c 

146 PR=0.711 

147 SC=0.595 

148 (e 

149 c 

150 c SEPARATE HEAT TRANSFER TERMS CALCULATED 
151 Cc 

152 c 

153 QC=(TS-TMEAN) /LHS 

154 Cc 

155 C 

156 c 

157 6 

158 IF(WET.EQ.0.) GO TO 50 

159 Cc 

160 ce 

161 Cc A POLYNOMIAL EXPRESSION (LOWE,1976) USED 
162 c TO COMPUTE SAT. VAP. PRESS. 
163 c : 

164 AO=6. 107799961 

165 A1=4.436518521E-01 

166 A2=1.428945805E-02 

167 A3=2.650648471E-04 

168 A4=3.031240396E-06 

169 c 

170 T=TMEAN 

7A EATA=AO+T *(A14+T*(A2+T*(A3+A4*T) ) ) 
shes C 

173 S0=676 .967 

174 S1=-0.0680427 

175 $2=-0. 120855E-02 

176 S$3=-0.518213E-05 

177 S$4=-0.249304E-07 

178 Cc 


179 TM=TS 
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LHW=SO+TM* (S1+TM*(S2+TM*(S3+S4*TM) ) ) 


LHW=(LHW/O. 2389) *1.0E07 
EPS=0.622 

QE 1=(PR/SC)**(0.63) 
QE2=-(EPS)/(PMEAN*CP ) *(LHW/LHS) 
QE3=EATA 

QE=QE1*QE2*QE3 


GOMTOr24 
QE=0.0 


CALCULATE REYNOLDS NUMBER FOR NEW DIAMETER. 


RE=TV*NDIAM/NU 


USING THE SUBLIMATION EXPRESSION DERIVED COMPUTE CHANGE 
IN MASS,HENCE CHANGE IN DIAMETER, AFTER TIME DELTA. 


NUSSELT RELATION ACCORDING TO ZUKAUSKAS 
APPROPRIATE FOR REYNOLDS RANGE 10**3-10**5 


NUSELT=0.23*(RE**0.6) 


CONST=(2*KA*NUSELT ) /RHOCO2 

DELD=( (CONST/NDIAM) *(QC+QE) ) *DELT 

ND IAM=NDIAM+DELD 
NMASS=PI*RHOCO2*CYLL*( (NDIAM*NDIAM) /4) 


IF(NDIAM.LE.O.) GO TO 25 
GO TO 30 
NDIAM=0.0 


WRITE(7,14) TIME,TA,PA,H,NDIAM 
IF(NDIAM.EQ.0.) GO TO 10 
IF(H.LE.3000.) GO TO 10 

GO TO 20 


STGP 
END 
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